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Lattice relaxations, surface phonon spectra, surface en- 
ergies, and work functions are calculated for Rh(lOO) and 
Rh(llO) surfaces using density-functional theory and the full- 
potential linearized augmented plane wave method. Both, 
the local-density approximation and the generalized gradi- 
ent approximation to the exchange-correlation functional are 
considered. The force constants are obtained from the di- 
rectly calculated atomic forces, and the temperature depen- 
dence of the surface relaxation is evaluated by minimizing the 
free energy of the system. The anharmonicity of the atomic 
vibrations is taken into account within the quasiharmonic ap- 
proximation. The importance of contributions from different 
phonons to the surface relaxation is analyzed. 
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I. INTRODUCTION 

Experimentally, the structure of a selid surface can 
be determined by LEEDH, ion scatteringj. X-ray diffrac- 
tionH, and helium-atom scatteringH. Surface phonons can 
be measured by high-resolution electron energy-loss spec.- 
troscopy (EELS)cl and inelastic helium-atom scatteringB. 
It is well understood that the first interlayer separation 
of metal surfaces is typically contracted at room temper- 
atureQ, and only few exceptions to this "rule" have been 
discovered so far. Due to the anharmonicity of the inter- 
atomic potentials, both the surface interlayer spacing and 
the bulk lattice constant depend on temperature. Fur- 
thermore, because of lower symmetry, the anharmonicity 
at the surface is expected to be more significant than that 
in the bulk. |-| ■— ■ 

First-principles calculations oOattice relaxationsETtll 
and surface phonons for metalal3^ll3-|typically employ 
the density-functional theory (DFT)EJ. In most previ- 
ous work, the equilibrium geometry was determined by 
minimizing the total energy of the system (we will call 
this the "static equilibrium"). We note that the surface 
structure obtained in this way corresponds to T = K 
and neglects the effects of zero-point vibrations. With re- 
spect to the latter we note that for example for the crys- 
tal bulk zero-point vibrations give rise to an increase of 
the lattice constant by about 0.2-0.5%, compared tp^t 
geometry determined by the total-energy minimunJlJl 
Starting from the "static equilibrium" geometry, surface 
phonons were calculated either directly (using the su- 
percell approach)t3TiiL-ar by applying density-functional 



At finite temperature, the equilibrium geometry is de- 
termined by the minimum of the free energy, which, in 
addition to the total energy, includes the contribution 
from lattice vibrations. Typically this gives rise to a lat- 
tice epinansion, but it is also possible that contractions 
occurta. Because anharmonicity is more pronounced at 
the surface than in the bulk an increased ittijetlayer spac- 
ing of the top layer, di2, is-±o be expectedEJO. Recently, 
Narasimhan and ScheffierEj and Cho and ScheffleicJ per- 
formed a theoretical study of the thermal expansion of 
Ag(lll) and Rh(lOO). The free energy was calculated by 
including the contrihutimi of phonons in the quasihar- 
monic approximationtjO. The temperature dependence 
of di2 was obtained by minimizing the free energy tbp 
considered slab with respect to di2- In these studiesotj, 
as a first approximation, only the top layer was allowed 
to vibrate as a whole parallel and perpendicular to the 
surface. The effective vibrational "modes" considered in 
these works correspond to modes at the Brillouin zone 
center for the top layer vibrating on a rigid substrate. No 
information about the true surface phonon spectrum was 
obtained, and thus, not even an attempt was made to dis- 
cuss how a good-quality summation over the phonon Bril- 
louin zone may affect the result. In fact a better treat- 
ment is quite elaborate as it requires the diagonalization 
of the dynamical matrix of the whole slab. In order to 
improve on the previous calculationscj, we present in this 
paper a theoretical study of the structure and dynamics 
of Rh(lOO) and Rh(llO) surfaces by including vibrations 
of the whole slab. The interplanar force constants are 
obtained from the directly calculated atomic forces using 
the fuU-potentia] linearized augmented plane wave (FP- 
LAPW) methodkaka. It has been shownlBlB that the 
surface phonon modes for wave vectors at high-symmetry 
points of the surface Brillouin zone (SBZ) can be deter- 
mined from the knowledge of interplanar force constants. 
In the present work, these force constants and the corre- 
sponding surface phonons are calculated as a function of 
the interlayer spacing at the surface. The temperature 
dependence of surface relaxation is then determined by 
minimizing the free energy within the quasiharmonic ap- 
proximation. Contributions from different phonon modes 
to the vibrational free energy will be examined. In Sec- 
tion |ll[ wc will describe some details of the calculation 
method. The results are presented in Section III. 
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II. METHOD OF CALCULATION 
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A. General theory 

We consider a crystal slab consisting of a finite number 
of atomic layers perpendicular to the z direction, and 
of infinite extension in the x and y directions. In the 
quasiharmonic approximation, the free-eujecgy function 
of the considered system can be written asliS 

F({dy},r) = <i>({d„}) + Fvib({d.,}-7^) 

where $ is the static total energy which can be obtained 
by first-principles calculations, fee and h are the Boltz- 
mann and the Planck constants, and {dij} is the set of 
inter-layer distances, di2, (i23, • ■ ■, between layers 1 and 
2, 2 and 3, etc. The vibrational free energy is denoted as 
Fvib, and ti^p(q||, {rfij}) is the frequency of the p-th mode 
for a given wave vector q|| , evaluated at the geometry 
defined by {dy }; and iV is the number of atoms in the 
slab. Anharmonicity of the inter-atomic potentials is in- 
cluded in this description because $ and the vibrational 
frequencies depend on the inter-layer distances {dij\. 
The free energy and the equilibrium distances {d^j} at 
a given temperature are determined by the minimum of 
F({dij},r) with respect to {dy}. 

In the present work, we have calculated the vibrational 
frequencies_and the corresponding free energy using the 
slab modelEIl and density-functional theory. We expand 
the static total energy, of the system in a Taylor series 

$ = $o + 5I'^«(l)""(l) 

l,a 

+ ^EE'^"/3(l'l>a(lW(l') + - , (2.2) 

\,a l',/3 

where Mq(1) is the a component (a — x,y,z) of the 
displacement of the 1-th atom from its mean position 
Ro(l) (taking the thermal expansion of the lattice into ac- 
count). The instantaneous position of an atom is given 
by r(l)=Ro(l)-l-u(l). We consider here the monatomic 
lattice, and the set of integers l=(Zi, ^2, ^a), which specify 
a particular atom, has the following meaning: I3 labels 
the crystal planes lying parallel to the surface, and h, I2 
specify the points in the two-dimensional lattice which 
spans a plane. In the quasiharmonic approximation the 
equations of motion are 

M^u^{\)^-Y,cp^p{\\')up{\') , (2.3) 

l',/9 

where M is the atomic mass, 0^^(1,1') is the force con- 
stant which is defined by 



The subscript "0" in Eq.(2.4) indicates that the force 
constants (ji^p (1, 1 ) are to be evaluated at the mean po- 
sitions of the atoms, rather than at the positions of the 
static equilibrium. Due to the two-dimensional transla- 
tional property of the slab, the normal mode solutions to 
Eq.(2.3) have the formEZl 

M„(l) = M-^'^Q^Co.{h)eMM\\ ■ Ro||(l||) 

+ q|l .Ro|,(Z3)-^i]} , (2.5) 

where Qq is the vibrational amplitude, Cq(^3) is a polar- 
ization vector which will turn out to be the eigenvectors 
of the dynamical matrix, and u is the-, vibrational fre- 
quency. Following the notation of Ref.EZi , we have q|| = 

{qx,qy), Row = (a;, ?/), Ro|| (1) = R-o||(l||) + R-olK's), 1|| = 
{hjh)- Insertion of Eq. (2.5) into Eq. (2.3) leads to the 
eigenvalue equation 

Y^D^p{kJ'3;qii)Ul'3;ci\\)^io\q\l)Cf3{h;q\l) , (2.6) 

where the elements of the dynamical matrix are defined 
by 

Da0{l3j3,<l\\) = J^^'PapikJ's^'^'w " 1||) 

xexp{iq|| • [Ro||(l'| - 1||) 
+ Ro||(/;,)-Ro||(?3)]} . (2.7) 

Once the force constants are obtained, the frequencies 
and polarization vectors for all modes in the slab for a 
given q|| can be obtained by diagonalization of a 3A^ x 3A^ 
matrix. 



B. Density functional theory calculation 

In the present paper, the considered systems Rh(lOO) 
and Rh(llO) are modeled by a periodic slab consisting 
of 7 layers of Rh and a vacuum region with the same 
thickness. We employ density-functional theory and 
the FP-LAPW methodOtJ. The exchange-correlation 
functicmal is treated by the local-density approximation 
(LDA)E3 aS|-aj"ell as the generalized gradient approxima- 
tion (GGA)cJ. We first calculate the lattice relaxations, 
surface energies and work functions of Rh surfaces at 
T = K neglecting the influence of zero-point vibrations. 
The geometry is optimized by a damped molecular dy- 
namics allowing the top two layers on both sides of slab 
to relax. The remaining atoms are kept at the bulk lat- 
tice sites. The energy cut-off for the FP-LAPW basis 
is taken to be 15 Ry and for the wave functions inside 
the muffin-tin spheres angular momenta are taken into 
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account up to = 8. The cut-oflF energy for the po- 
tential is taken as — 100 Ry. For the k||-sanipHng, we 
use a uniform mesh of 15 points in the irreducible part 
of SBZ of the (1 X 1) surface displaced from F. 

In order to calculate phonon frequencies Q£-.the slab, 
an extension of the "frozen phonon" methodllj is used. 
For a selected q|| in the surface Brillouin zone, we can 
calculate all the interplanar force constant matrices by 
distorting the slab in an appropriate way. The obtained 
force constants allow us to set up the dynamical matrix 
and solve for the eigenvalues and eigenvectors of the sys- 
tem. The interplanar force constants coupling planes Z3 
and ^3 are given by 



(fi'ap{h^'3.^,^.\\) = y^0a/3(^3,^3il'| - 1||) 



In 



exp [iqii • Ro||(l|| - 1||)] 



(2., 



In the harmonic or quasiharmonic approximation, the 
interplanar force constants, 0^^(^3, ^3; qy), are related to 
the forces by 



9fa(/3,q||) Af„G3,q||) 

V^afll'S, '3! ^11 j — Tj TTT ~ Tj JTT- 



, (2.9) 



where AFQ(/3,q||) is the a component of the force dif- 
ference at atomic layer ^3 under a distortion of layer 
in accordance with q||, m^(1||,Z3) is the /3 component of 
atomic displacement at (1||,Z3). The displacement for 
other atoms in layer Z3 is given by 



u(lj|,;^) = u(l||,;^)exp[iq||.Ro||(l||-l||)] 



(2.10) 



The dynamical matrix of Eq.(2.7) can then be written as 

^a/3(^3,^3;q||) = ]^€/3(^3,?3;q||) 

exp{zq|| • [Roll (/^,)-Ro| 1(^3)]} ■ (2.11) 

In order to calculate the thermal expansion, we need 
to know the interplanar force constants as a function of 
the interlayer distance. In most cases, the thermal ex- 
pansion of the lattice at the surface along z direction is 
expected to be more significant than the expansion of the 
bulk lattice constant. The latter also defines the expan- 
sion of interatomic distance in the x-y plane. We there- 
fore calculate the surface interplanar force constants as 
a function of di2. Furthermore, we calculate the inter- 
planar force constants for the crystal bulk. From all the 
interplanar force constants, the dynamical matrix (as a 
function of di2) of a thick slab can be constructed. In 
fact, for the phonon calculations it is important to use 
a significantly thicker slab than for the electronic struc- 
ture calculations, in order to obtain a reliable descrip- 
tion of surface modes, surface resonances and the surface 
phonon density of statesE3tZl. We therefore calculate the 
surface phonon spectra for Rh(lOO) and Rh(llO) for a 
201-laycr thick slab. The surface relaxation di2{T) is 
then obtained by minimizing the free energy of this slab 
for various temperatures. 



III. RESULTS 

For Rh bulk the theoretically obtained lattice con- 
stant (neglecting the influence of zero-point vibrations) is 
a^^ = 3.79 A for DFT-LDA and a^^ = 3.83 A for DFT- 
GGA. The FP-LAPW parameters employed for these 
calculations are the same as those for the surface cal- 



culations listed in Section II B, except that for the bulk 



study we use a mesh of 72 k-points in the irreducible 
part of the fcc_Brillouin zone. The experimental value is 
^cxp _ 3.80AE3. The surface properties of Rh are calcu- 
lated with the theoretical bulk lattice constant. In Table 
^ we present the results of surface relaxations, work func- 
tions and surface energies of Rh(lOO) and Rh(llO) cal- 
culated in the LDA and GGA neglecting the zero-point 
vibration. Results from other calculations and the ex- 
periments are also listed for comparison. The present 
results are in gjSjpd|-ajppement with those of the previ- 
ous calculationsoliiliHdS. In the following calculations of 
the surface dynamics, we will use the GGA. The surface 
interplanar force constant will be calculated as a func- 
tion of di2, while the interlayer spacing ^23 is kept as the 
static equilibrium value. 

The calculated interplanar force constants coupling the 
surface layer with other layers in Rh(lOO) slab are given 
in Table O. Three high symmetry surface wave vectors 
are considered (see Fig. |l|): F, X and M. These results 
are obtained by distorting the slab around the zero-force 
geometry at T = K, which is given in Table |. The cor- 
responding results for the interior of the slab are given 
as well. The interplanar force constant matrices for the 
Rh(llO) surface are given in Table [II. Four wave vectors 



at F, X, Y, and S points in the SBZ are considered (see 
Fig. 1^) . It is not surprising to see in Table ^ and Table HI 
that the intralayer part of the force constants ((/"^^(l, 1)) 
are most significantly changed at the surface compared 
to that in the bulk because of the truncation of surface. 
At the same time, the interlayer part of force constants 
at surface (^^^(1, Z3), ^3 = 2,3) are also modified due to 
the surface relaxations. For most of the matrix elements 



in Table III , there is a trend for an enhancement in mag- 
nitude for force constants coupling the surface layer to 
its neighbors compared with the bulk force-constant ma- 
trix elements. This is in agreement with the findings of 
Al(110)y. 

Anharmonicity is included in our studies by calculat- 
ing the force constants as a function of the interlayer 
distance di2. Figure || shows the variation of the surface 
interplanar force constants at the X point for Rh(lOO). It 
can be clearly seen that the interlayer force constants de- 
crease monotonically as di2 changes from 9% contraction 
to 3% expansion. This "softening" of the force constants 
reflects the anharmonicity of the bond strength at this 
surface. We also note that the intralayer force constants 
of the second layer ((/)^^(2, 2) ) are also significantly soft- 
ened at the X and M points for Rh(lOO), and at the X, Y, 
and S points for Rh(llO), with the increase of di2. While 
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at r, there are no intraplanar vibrations, the "softening" 
of 0^^(2,2) is already included in the surface interpla- 
nar force constant (/)^^(1,2). In the present work, the 
changes of the intraplanar force constants for the second 
layer with di2 are taken into account. 

Having obtained the force constants for different di2, 
we calculate the phonon frequencies of the slab for each 
di2. Then, using Eq. (2.1), the temperature dependent 
top-layer relaxations for Rh(lOO) and Rh(llO) are de- 
termined by minimizing the free energy with respect to 
di2. Figure m shows the variation of Adi2/do for Rh(lOO) 
with temperature in the case in which vibrations are cal- 
culated from different wave vjactors qy in the SBZ. The 
results of Cho and Schefflem are also given for com- 
parison. It can be seen that when only the F point vi- 
bration is taken into account, the variation of Adi2/do 
with teiaperature is in close agreement with the results 
of Ref.Ej. The small difference is due to the fact that 
only the top layer was allowed to vibrate in Ref.E3, while 
in the present work, vibrations of the whole 201-layers 
thick slab are included. At low temperature (below 300 
K), the vibrations from the F, X, and M points contribute 
with similar importance to the surface expansion. The 
differences between the contributions of the different qy 
vectors becomes significant at high temperature because 
of the different frequency distributions between F, X and 
M. 

Using the "frozen" phonon mcthodElil'0, it is not 
practicable to get the phonon frequency for an arbitrary 
qy point in SBZ. For phonon, in contrast to electrons, it 
is in fact acceptable to approximate the summation over 
the whole SBZ by the contributions from only the high- 
symmetry q|| points. The weights of F, X, and M are 
J, J, and i, respectively. Thus, the vibrational free en- 
ergy is given by F^ib = ^Fvib(r) + |i^vib(X) + iFvib(M). 
The temperature dependence of the top-layer relaxation 
for Rh(lOO) by including all the vibrational contributions 
from X, M, and F is shown in Fig. ^ . The arrow marks 
the result of the surface relaxation obtained when lat- 
tice vibrations (also zero-point effects) are neglected. We 
find that at low temperature the tliermal expansion for 
Rh(lOO) is close to the resuft of Ref.El. At 3QQK, the top- 
layer relaxation is —1.64%, the value of Ref.Ej is —1.45%. 
However, for higher temperature the thermal expansion 
will be overestimated if only the vibrations at F point are 
included. In general, aux improved calculations confirm 
the conclusion of Ref.tJ. The calculated surface ther- 
mal expansion coefficient is as — [di2)~^ {dd\2 / dT) = 
40.7 X lO-^K-i at T=300 K, which is 5^0. times larger 
than in the bulk (a^ = 8.2 x IQ-^K-iB). Thus for 
Rh(lOO) the thermal expansion is clearly increased at the 
surface. The different contributions to the surface relax- 
ation from vibrations parallel to the surface and perpen- 
dicular to the surface are also shown in Fig. IJ. All the 
vibrational modes at X, F, and M are included. The re- 
sults in Fig. ^ confirm that the parallel vibrations have 
a significantly larger influence on the surface thermal ex- 



pansion than the perpendicular vibrations (z-direction). 
This g g agreement with the findings and conclusions of 

We now turn to the more open (110) surface. The 
temperature dependence of the top-layer relaxation for 
Rh(llO) is shown in Fig. |5[ The vibrational free energy 
is calculated by summing over the wave vectors F, X, Y, 
and S in the SBZ. These four points all have the same 
weight, therefore the vibrational free energy is given by 
i^vib = i[i^vib(r) + Fvib(X) + i^vib(Y) + f;ib(S)]. Similar 
to Rh(lGG), the contribution to the thermal expansion 
of c?i2 from in-plane vibrations is larger than that from 
out-plane vibrations. At r=300 K we obtain the sur- 
face thermal expansion coefficient ag = 59.4 x IG^^K^^, 
which is 7.2 times larger than the bulk thermal expansion 
coefficient. Comparing the thermal expansion coefficient 
of Rh(llO) with that of Rh(lOO), we find that Rh(llO) ex- 
hibits more significant anharmonicity than Rh(lOO). The 
top-layer relaxation is Adi2/do — —7.4% at r=300 K 
which is_dose to the result of a room temperature FEED 
analysia^, which gave a value of Afii2/fio = —6.9%. The 
surface relaxation obtained when just the total energy 
is minimized is noticeably different, namely Adi2/c?o = 
-9.2%. 



In Table LV, we present the calculated surface phonon 



frequencies corresponding to the surface relaxation at 
T=300 K. The local density of states are obtained from 



^q||,a(w) = ^(5(w-Wp)|Ca(l,q||,p)|^ 



(3.1) 



where Ca{^,(\\\,p) is the surface layer polarization vector 
for the p-th normal mode of the slab with respect to wave 
vector q||, a denotes x^y and z. The calculated surface 
phonon densities of states are shown in Fig. || and Fig. |^. 
Different plots correspond to different polarizations. At 
the M point for Rh(lOO) ( Fig. results for y polariza- 
tion arc identical to those for the x polarization due to 
symmetry. It can be seen that in all cases (Fig. ^ and 
Fig. 0), the surface modes are strongly localized at sur- 
face and are well separated from the bulk modes. Thus it 
should be possible to identify them in EEFS or He scat- 
tering experiments. It is these localized surface phonons 
that are most sensitive to the surface relaxations, and 
hence govern the anharmonic thermal expansion at the 
surface. It can be seen in Figs. || and ^ that there are 
more localized surface phonon modes polarized in the x 
and y directions than in the z direction. This is one of the 
reasons why the parallel vibrations give a more important 
contribution to the thermal expansion (see Figs. ^and|^). 



IV. SUMMARY 

In conclusion, we have used full-potential FP-FAPW 
method to investigate the surface structures and surface 
phonon spectra of Rh(lOO) and Rh(llO) surfaces. We 
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studied the thermal expansion by taking anharmonic- 
ity into account within the quasiharmonic approxima- 
tion. It is found that different phonons give differ- 
ent contributions to the surface relaxation in the range 
of high temperature (higher than 300 K). The surface 
thermal expansion coefficients calculated by including 
all the high symmetry points in surface Brillouin zone 
are 40.7xlO-<^K-i for Rh(lOO) and 59.4xlO-*'K-i for 
Rh(llO) at T=300 K, which are 5.0 and 7.2 times larger 
than that of the bulk, respectively. The obtained surface 
relaxations of Rh(lOO) and Rh(llO) at T = 300 K are in 
agreement with the experimental measurements at room 
temperature. The calculated results confirm that the in- 
planc vibrations give a more important contribution to 
the surface thermal expansion than the out-of-plane vi- 
brations. 
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Figure 1 




FIG. 1. Surface Brillouin zones for fee (100) and (110) sur- 
faces. High symmetry points are indicated. 
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FIG. 3. Top-layer relaxation of Rh(lOO) as a function of 
temperature calculated in different approximations, i.e., in- 
cluding different vibration modes. The results of Cho and 
SchefHercJ are shown for comparison. 
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Figure 2 

FIG. 2. Variation of the surface interplanar force constants 
at X with Adi2/do for Rh(lOO). 



6 



T3 

< 



Rh(100) 



result of Cho et al. 
only perpendicular modes 
only parallel modes 
all modes 




-5.0 



100 200 300 400 
Temperature (K) 



500 600 



-10.0 



Rh(110) 

A all modes 
O only parallel modes 
• only perpendicular modes 




100 200 300 400 500 600 

Temperature (K) 



Figure 4 



Figure 5 



FIG. 4. Top-layer relaxation of Rli(lOO) as a function of 
temperature calculated by including the vibration modes at 
r, X, and M. The contributions from vibration parallel and 
perpendicular to the surface are are displayed. The arrow 
points at the result obtained when vibrational contributions 
to the free energji-|are neglected completely. The results of 
Cho and ScheffieiO are shown for comparison. 



FIG. 5. Top-layer relaxation of Rh(llO) as a function of 
temperature calculated by including the vibrational modes at 
F, X, Y, and S. The contributions from vibrations parallel 
and perpendicular to the surface are displayed. The arrow 
points at the result obtained when vibrational contributions 
to the free energy are neglected completely. 
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Figure 6 



FIG. 6. Surface phonon density of states of Rh(lOO) for 
the wave vectors at the X and M points in the surface Brillouin 
zone. 
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TABLE II. Interplanar force constants of Rh(lOO) cou- 
pling the surface layer to other layers for wave vectors at F, 
X, and M at Adia/do = -3.0%, Adas/do = -0.1%. The 
units are 10* dyn/cm. The corresponding interplanar force 
constants for the interior layers of the slab are also shown for 
comparison. The indices (a/3) and the arguments (hjl's) of 
(I>^0{l3, 1's) are combined as {hajl'^P). Matrix elements not 
listed are zero by symmetry. 

{ha, I'^P) Surface Interior {ha, I'-^fd) Surface Interior 









qii =r 






{lx,2x) 


-4.89 


-4.92 


{Ix, 3x) 


-0.03 


-0.30 


{lz,2z) 
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-7.65 


{lz,3z) 


-1.45 


-1.84 
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26.23 
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0.53 
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10.18 


(iy,3j/) 
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0.03 


{Iz, Iz) 


10.15 


17.07 


{lz,3z) 


-0.57 


-0.70 


{lx,2z) 


-6.42 


-5.46 


{Iz, 2x) 


-5.41 


-5.46 
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Figure 7 

FIG. 7. Surface phonon density of states of Rh(llO) for 
the wave vectors at X , Y, and S points in the surface Brillouin 
zone. 



TABLE I. Surface relaxations Adia/do and Adas/do, work 
function (j) (eV) and surface energies a (eV/atom) for Rh(lOO) 
and Rh(llO) surfaces obtained by different calculations and 
experiments. The results correspond to T = K, and none of 
the theoretical values includes the effect of zero-point vibra- 
tions. 



Rh(lOO) LD 

LDAlii 
LDAp 



Adi2/do 



Adaa/do 



ggaE 

this - LDA 
this - GGA 



-3.5% 
-3.8% 
-3.0% 
-2.8% 
-3.3% 
-3.0% 



experiment -|-0.5±1.0%Ji 

experiment — 1.16±1.6%E 

Rh(llO) LDA2£ -10.2% 

LDAliil -9.8% 

this - LDA -9.9% 

this - GGA -9.2% p 

experiment -6.9±1.0%fi 



-fO.7% 
-0.2% 
-0.1% 
-0.2% 
-0.1%^ 
0±1.5%^ 
0±1.6%E^ 

4-2.5% 
-f2.6% 
-f2.0% 
+2.l7o r 
-hl.9±1.0%E 



5.25 

5.26 
4.92 
5.30 
4.9, 
4.6, 



1.27 
1.44 
1.29 
1.04 
1.23 



TABLE III. Interplanar force constants for Rh(llO) at F, 
X, Y, and S at Adia/do = -9.2%, Adas/do = 2.1%. The 
units are 10* dyn/cm. The corresponding interplanar force 
constants for the interior layers of the slab are also shown for 
comparison. The indices (a/3) and the arguments (/3,/3) of 
4>^^l){h,l'-i) are combined as {ha,l'-j/3). Matrix elements not 
listed are zero by symmetry. 



3i i.iM 
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qii = Y 






{Ix, Ix) 


5.69 
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TABLE IV. Calculated surface phonon frequencies v for 
Rh(lOO) and Rh(llO) corresponding to surface relaxation at 
T=300 K. Units arc in THz. 
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